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Abstract 

Relationships between the coupling constant and the binding energy of threshold bound states 
are obtained in a simple manner from an iterative algorithm for solving the eigenvalue problem. 
The absence of threshold bound states in higher dimensions can be easily understood. 
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In a mathematically elegant paper, Simon studied the one and two dimensional 
Schrodinger operators —d 2 /dx 2 + XV (x) and — A +AV(x) where either V(x) or V(x) is 
described in this work as the potential and the parameter A (> 0) is termed the strength 
of the potential. Simon provided necessary and sufficient conditions for the existence of a 
bound state when A is small In one dimension a threshold bound state (i.e., one just 



bound) exists 
function of A|2 



'or many finite short-range potentials and its binding energy is an analytical 



[J. Furthermore, using the theory of trace class determinants 0, Q|, a simple 
expansion for the binding energy of the threshold bound state has been obtained (see Ref[lfl). 
More recently, Gat and Rosenstein have pointed out that perturbative methods provide a 
suitable means for calculating the binding energy of this state jsl]. This is somewhat peculiar 
since a convergent expansion for the binding energy in A exists, but no apparent poles appear 
in the expansion of the S matrix to any finite order in perturbation theory. Rather than 
use perturbation theory, we wish to point out that the expression for the binding energy 
of the threshold bound state obtained by Simon in one and two dimensions can easily be 

n 

obtained from a simple non-perturbative iterative algorithm[7| and we provide an intuitive 
explanation of the form of the expansion in different dimensions. 

In the algorithmic approach, eigenvalues and the associated eigenfunctions are determined 
as functions of the strength of the potential, A. To illustrate the method, we consider the 

n 

one-dimensional eigenvalue equation [7| 

[-d 2 x - XV{x)]u{x) = -eu(x) (1) 

subject to 

lim u{x) = 0. (2) 

\x\ — >oo 

Here d x = d/dx, A > and f V(x)dx > 0. We shall always assume V(x) — > as |x| — > oo. 
The energy eigenvalue, — e (with e > 0), is negative and corresponds to a bound state. Using 
Green's method a solution to equation (pQ) is given by 

/oo 
G e (x — x')V{x')u{x')dx' (3) 
■oo 

where the Green's function G e (x) satisfies 

[-dl + e]G e (x)=5(x) (4) 
lim G e (x) = 0. (5) 

\x\ — >oo 



Normalizing u(x) at an arbitrary x re f 

u(Xref) = 1 (6) 

allows A to be written as (see Eq. Q) 

A = f™ oo G £ (x ref -x')V(x')u(x')dx' (7) 

which can then be used to eliminate A from equation (J3j) : 

f°° GJx - x')V(x')u(x')dx' 
G e (x re f — x')V(x')u(x')dx' 

Using equations ([7]) and @, A can then be determined as a function of e as follows. For a 
particular choice of e Eq. ((8]) can be iterated from a reasonable starting point, uq(x): 

G e {x - x')V{x')u n {x')dx' 
n+1 J^ao G e (x ref - x')V(x')u n (x')dx' 

until it converges and A can then be determined from Eq. ([7j). Repeating this procedure for 
a different value of e yields a different value of the strength of the potential, A. When enough 
(e, A) pairs have been determined, a simple interpolation procedure can be used to determine 
the dependence of e on A. Furthermore, for larger values of e a simple relationship between 
A and e can be obtained for non-singular symmetric potentials which vanish asymptotically 
which can be used to make the the algorithm more efficient ^|. 

For small values of e, corresponding to states on the threshold of being bound, the an- 
alytical dependence of e on A may be obtained approximating the Green's function, which 
satisfies Eq. {4]) and is given by 

e ~Ve\x\ 

o.(«) = ^r- do) 

Expanding G e (x) in e: 

G e (x) = ^= + ... (11) 
and substituting this into Eq. ([Sj) yields 

u(x) = l + .... (12) 

to leading order in e. From Eq. (£3) one therefore easily obtains the following approximate 
relationship between the coupling constant A and e 



which is valid for small values of e. For arbitrary small coupling it provides an analytical 
expression for the lowest bound state for a large class of potentials in one dimension provided 
J' x oo V(x)dx > o[lj]. Furthermore, had we included the standard factor of 1/2 in the first 
term of the eigenvalue equation, Eq. (00), then we would obtain 

i / ^oo \ 2 

1 \2 



2 X'U V(x)dx) (14) 

which is precisely the result obtained by Simon to 0(A 2 ). 

An important feature of Eq. ([3]) is that it has a solution for e, for arbitrarily small A. 
This follows since the Green's function, Eq. (fTUI) . is unbounded from above as e — ► 0: 

lim G € (x) = +oo. (15) 

Hence for very small values of A, the value of e can always be adjusted until the product 
XG t (x) is non-negligible. 

We note these results can be generalized to higher dimensions since it can easily be seen 
that Eqs j3|) and (JSj) become 

u(x) = A I G e (x - x')^(x>(x')cfV (16) 

J all svace 



and 



, , l, space q(x-x^(x>(xVV 

1 ) Lls paC eGe{*r e f-x!)V{*)u{*)d^- 1 } 

The essential difference arises only from the different form the Green's function takes in 
different dimensions. 
In two dimensions 

G e (x) = ^o( e |x|) (18) 

a 

where K Q (-) is a Bessel function of the second kind of order zero[9||. Expanding _K" (e|x|) for 
small e yields 

#o(e|x|) = In (1/e) + In {2e^/\x\) + 0(e 2 ) (19) 

where 7 = 0.57721... is Euler's constant. Thus from Eq. (fPTj) we obtain, for sufficiently 
small e, 

u(x) = 1 + . . . (20) 

and therefore 

^ — 7T\ r~7 (21) 



As in one dimension provided Jv(x)d 2 x > 0, a threshold bound state exists for arbitrarily 

small A if J \V(x)\ 1+l3 d 2 x < oo (some (3 > 0) and J (1 + x 2 f\V{^)\d 2 x < oofl. Again it is 
the divergence of the Green's function, at fixed spatial argument, when e — > 0, that leads to 
a threshold bound state at arbitrarily small A. 

On the other hand in three dimensions the Green's function is 

G ' (x > = <22) 

At fixed |x|, this does not diverge as e — > 0. This is suggestive of the known fact that 
in three and higher dimensions an arbitrarily weak attractive potential does not possess a 
bound state jl]; there has to be a certain strength of the potential before it can support a 
bound state. We note that the leading term in an expansion of e of the Green's function, in 
n = 3 and higher dimensions, is not independent of x. This is different to the corresponding 
behaviour of the Green's functions when n = 1 and n = 2 and suggests that n — 1 and 
n = 2 which may be thought of as being atypical of all other dimensions. 

To understand the property of an arbitrarily weak attractive potential to bind a particle 
in n = 1 and n = 2 two dimensions but not in three or more dimensions, we can relate it to 
an apparently different problem of how much time a random walk in n dimensions spends in 
the vicinity of its starting position. Using Dirac notation in the general n-dimensional case 
where p is the momentum operator and |x) ((x|) is an eigenket (eigenbra) of the coordinate 
operator, we have 

roc 

G e (x) = (x|(p 2 + e)- 1 ^) = / dt(x|e-(P 2+£ )*|0). (23) 

Jo 

This quickly leads to 

roc -|x| 2 /(4t) 

limG e (x)= / T- , w dt. (24) 
e^o el ; J (Ant)-/ 2 K ' 

In a formulation of the random walk in discrete space and discrete time[l0|, it is precisely 
the analogue of lim e ^ G e (x) which determines the mean time a random walk spends in the 
vicinity of a site at position x, given it was at position at time t — 0. In n = 1 and n = 2 
dimensions lim^oG^x) is infinite, implying an infinite amount of time is spent at x. By 
contrast, when n > 3, Go( x ) is finite. In the continuous space quantum mechanical problem 
considered here, it is precisely the finiteness (or lack of finiteness) of lim e ^ G e (x) that 
determines the dimensionalities where an arbitrarily weak potential can possess a threshold 
bound state. 
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We conclude by pointing out that at least for small values of A, good approximate ana- 
lytical relationships between e and A exist in one and two dimensions which may be used to 
improve the convergence rate of the aforementioned iterative algorithm 7]. 



[1] B. Simon, Ann. Phys. (N. Y.) 97, 279 (1976). 

[2] S. Fliigge, Practical Quantum Mechanics, vol. 1 (Springer- Verlag, Berlin, 1970). 
[3] M. Znojil, J. Phys. A 20, 1751 (1987). 

[4] N. Dunford and J. Schwartz, Linear Operators, Part II: Spectral Theory (Interscience, New 
York, 1963). 

[5] I. C. Gohberg and M. G. Krein, Transl. Amer. Math. Soc. 18 (1969). 

[6] G. Gat and B. Rosenstein, Phys. Rev. Lett. 70, 5 (1993). 

[7] D. Waxman, J. Phys. A 31, 1329 (1998). 

[8] W. A. Berger and H. G. Miller, J. Phys. A 39, 14075 (2006). 

[9] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions (Dover, New York, 1965). 
[10] C. Itzykson and J. Drouffe, Statistical Field Theory, vol. 1 (Cambridge University Press, Cam- 
bridge, 1989). 



6 



